Fermionic massive modes along cosmic strings 
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The influence on cosmic string dynamics of fermionic massive bound states propagating in the 
vortex, and getting their mass only from couphng to the string forming Higgs field, is studied. Such 
massive fermionic currents are numerically found to exist for a wide range of model parameters and 
seen to modify drastically the usual string dynamics coming from the zero mode currents alone. In 
particular, by means of a quantization procedure, a new equation of state describing cosmic strings 
with any kind of fermionic current, massive or massless, is derived and found to involve, at least, 
one state parameter per trapped fermion species. This equation of state exhibits transitions from 
subsonic to supersonic regimes while the massive modes are filled. 



I. INTRODUCTION 



Since it was realized that some early universe phase transitions might lead to the formation of topological defects [0 , 
cosmic strings have been the subject of intense work within the context of cosmology [Q. The large scale structure 
, generated by an ordinary string network in an expanding universe, as well as its imprint on the microwave background, 
^ ' have thus been derived pLE [in order to state on their significance in the wide range of mechanisms in which they had 
^\ I been originally involved [^^. These predictions, compared with the observations therefore constrain the symmetry 
' breaking schemes effectively realized in the early Universe. These, associated with the most recent data for the 
microwave background anisotropics [Q, even seem to show that such ordinary string networks could not have play 
^5 the dominant role in the Universe evolution, thereby all the more so constraining the particle physics symmetries 
leading to their formation. However, as was recently shown a non-negligible fraction of such defects could have 
contributed to the overall cosmic microwave background (CMB) anisotropics. 

Meanwhile, it was shown by Witten |8j that in realistic physical models, involving various couplings of the string 
forming Higgs field to other scalar or fermion fields, currents could build along the strings, turning them into "super- 
conducting wires." Without even introducing couplings with the electromagnetic fields Q, the breaking of Lorentz 
O i' invariance along the vortex induced by such currents may drastically modify the string properties, and thus, the 
O , cosmological evolution of the associated networks. In particular, cosmic string loops can reach centrifugally supported 
1^ ' equilibrium state, called vortons that would completely dominate the Universe Theories predicting sta- 

ble vortons thus turn out to be incompatible with observational cosmology, hence the particular interest focused on 
"superconducting" models. 

^ . , Unfortunately, all the new properties and cosmological consequences stemming from string conductivity have not 
■ yet been clearly established, because of the complicated, and somehow arbitrary, microphysics possible in these 
models. However, although the microscopic properties induced by such currents depend on the explicit underlying 
field theory a macroscopic formalism was introduced by Carter Q which permits a unified description of 

the string dynamics through the knowledge of its energy per unit length U and tension T. These ones end up being 
functions of a so-called state parameter ui, as the current itself, through an equation of state. Such a formalism is, in 



particular, well designed for scalar currents, as shown in, e.g., Refs. |15 lq|: due to their bosonic nature, all trapped 



scalar particles go into the lowest accessible state, and thus can be described through the classical values taken by 



the relevant scalar fields |17|. The induced gravitational field |1J,|19[ or the back reaction effects |2C[| depend only 
on this state parameter. The classical string stability ||2l|,^ has already been investigated for various equations of 
state relating U and T, on the basis of scalar and chiral currents microphysics |^,^. Moreover, it was also shown, 
through a semiclassical approach, that fermionic current carrying cosmic strings, even though in principle involving 
more than one state parameter can also be described by an equation of state of the so-called "fixed trace" kind, 
i.e., U + T = 2M^. Such a relationship has the property of allowing stable loop configurations to exist, at least at the 
classical level |2^. Nevertheless, these results have been derived for fermionic currents flowing along the string in the 
form of zero modes only, as they were originally introduced by Witten , although it was shown that the fermions 
may also be trapped in the vortex with nonvanishing masses p6[ | : hence the following work in which the influence of 
such massive modes is studied for the simplest of all fermionic Witten model. 

In this paper, after deriving numerically the relevant properties of the trapped massive wave solutions of the Dirac 
equation in the vortex, we show that the quantization procedure, originally performed to deal with the fermionic zero 
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modes 25 1, can be generalized to include the massive ones, and leads to a new equation of state with more than one 



state parameter. In particular, it is found that the fixed trace equation of state, that holds for massless fermionic 
currents alone, is no longer verified. Besides, the massive modes are actually found to rapidly dominate the string 
dynamics, thereby modifying the classical vorton stability induced by the zero modes alone. 

Let us sketch the lines along which this work is made. In Sec. ^ the model and the notations are set, while we 
derive the equations of motion. Then, in Sec. Ill, by means of a separation between transverse and longitudinal 
degrees of freedom of the spinor fields, the massive wave solutions along the string are computed numerically for a 
wide range of fermion charges and coupling constants. The constraint of transverse normalizability is found to be 
satisfied only for particular values of the trapped modes mass, TtT say, whose dependence with the model parameters 
is investigated. The two-dimensional quantization of the m normalizable massive modes is then performed in Sec. IV, 
using the canonical procedure. In the way previously discussed in the case of zero modes |25| , the conserved quantities, 
i.e., energy-momentum tensor and charge currents, are then expressed in their quantum form. Their average values, 
in the zero-temperature case, and infinite string limit, lead to macroscopic expressions for the energy per unit length 
U and tension T which end up being functions of the number densities of fermions propagating along the string. Their 
derivation and extension to any kind and number of fermionic carriers is performed in Sec. [V|, while the cosmological 
consequences of this new analysis are briefly discussed in the concluding section. 



II. MODEL 

We shall consider here an Abelian Higgs model with scalar $ and gauge field B^, coupled, following Witten to 
two spinor fields, 'J and X say. Since we are only interested in the purely dynamical effects the current may induce 
on the strings, we will not consider any additional electromagneticlike coupling of the fermion fields to an extra gauge 
field. Thus, we consider here the so-called "neutral limit" ||l5|l 



A. Microscopic Lagrangian 

The previous assumptions imply one needs one local U{1) symmetry which is spontaneously broken through the 
Higgs mechanism, yielding vortices formation. The Higgs field is chosen as complex scalar field with conserved charge 
qcif, under the local U{1) symmetry, associated with a gauge vector field B^. The two spinor fields acquire masses from 
a chiral coupling to the Higgs field, and have opposite electromagnetic charges in order for the full (four-dimensional) 
model to be anomaly free ||^. Under the broken symmetry they also have conserved charges qc^^^, qcip-^ and qc-^^, 
qc^^^ for their right- and left-handed parts, respectively. With C^, Cg and C^, the Lagrangian in the Higgs, gauge, 
and fermionic sectors, respectively, the theory reads 

C = Ci, + Cg + C^+£x, (1) 

with 

^h = ^(i?;.<i>)t(l?^$)-m, (2) 



2'^ 
1 

4^ 



Cg^--H,,H^'^, (3) 



C^^'- [Xj^D^X - (D;j^)rx] - gX^-^X<S>* - gX^-^X<f, (5) 



where the U{1) field strength tensor and the scalar potential are 



H^, - V^B, - V,B^, (6) 

v{<i>) = ^m'-v'r, (7) 



while covariant derivatives involve the field charges through 
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and the relation 



D„X = V 



iq 



iq 



+ iq 



should hold in order for the Yukawa terms in £^ and to be gauge invariant. 



B. Basic equations 



(8) 
(9) 

(10) 



(11) 



This theory admits vortex solutions which are expected to form in the early universe by means of the Kibble 
mechanism jll. A cosmic string configuration can be chosen to lie along the z axis, and we will use Nielsen-Olesen 
solutions of the field equations p7|]. In cylindrical coordinates, the string forming Higgs and gauge fields thus read 



B, = B{r)5_ 



tie, 



(12) 



where the winding number n is an integer, in order for the Higgs field to be single valued under rotation around the 
string. In such vortex background, the equations of motion in the fermionic sector, for both spinor fields T read (here 
and in various places throughout this paper, we shall denote by T an arbitrary fermion, namely a spinor or X) 



ir^^J^^^B^ + M^J^ 



with the fermionic gauge currents 



3t = 1 



(13) 



(14) 



and the mass terms 



= gip cos nd + igip^5 sin n9, 



- gif cos — tgfj5 sm nb 



(15) 
(16) 



Note the fermionic currents have an axial and vectorial component because of the chiral coupling of the fermions to 
the Higgs field, as can be seen through the mass terms Mjr in Eqs. (15) and (p^). Moreover, since the Higgs field 
vanishes in the string core while taking nonzero vacuum expectation value, rj say, outside, the mass term acts as an 
attractive potential. As a result, fermionic bound states, with energy between zero and gr]^ are expected to exist and 
propagate in the string core. 



III. FERMIONIC BOUND STATES 



A. Trapped wave solutions 



Since the string is assumed axially symmetric, it is convenient to look for trapped solutions of the fermionic equations 
of motion, by separating longitudinal and transverse dependencies of the spinor fields. Using the same notations as 
in Ref. p^, the two-dimensional plane- wave solutions along the string, for both fermions, read 



( Ci(^)e"""^^ \ 



^2(r)e"""'^ 
6(^)e"""-'^ 



J 



-iliO 



fCiir 
V C4(r)e 



\ 



(17) 
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where e = ±1 labels the positive and negative energy solutions. Similarly to the Higgs field case, the winding numbers 
of the fermions, rrii and U, are necessary integers. In order to simplify the notations, it is more convenient to work 
with dimensionless scaled fields and coordinates. With mh = rf^/X the mass of the Higgs boson, we can write 



if = rjH, Q = n + qcd, B, and r = 

TOh 

In the same way, the spinorial components of the "if field are rescaled as 



(18) 



Voj + kai{g), £,2{g) = - ka2{g), 

271 vStt 



TOh 



\/uj - kas^g), ^4{g) = z^^Vw + ka4{g). 



(19) 



In the chiral representation, and with the metric signature (+,—,—,—), in terms of these new variables, Eqs. ( [T^ ) 
and dH) yield, for the ^ field. 



-i(mi — 1)6 



-i(m2 + l)e 



-i(m3-l)e 



-i(mi + l)e 



^ - .fi{Q)ai{g) 
ag 

^ - I2{g)a2{g) 

dg 



doa 
d^i 

da4 
dg 



f3{0)a3ig) 
fi{g)ai{g) 



s—e-'-^'0a2{g) - — ff(f?)e-*('"^+")^54(f^), 

TOh TOh 

_^^g-^„u9g^(^) + !!^^(^)e-(™3H-")«53(g), 

TOh TOh 

_^^e-™^«54(^^) + — H(f?)e'^("^-")^52(f^), 

TOh TOh 



(20) 



m 



= e- 



nif 



rrih 



^3(9) H{g)e 



-i{ra-i—n)6 



TOh 



ai{g). 



where mf grj is the fermion mass in the vacuum in which the Higgs field takes its vacuum expectation value 77, 
and m = \/ uP- — k'^ is the mass of the trapped mode. The coupling to the gauge field appears through the purely 
radial functions /: 



7 , s c^^Q-n mi ~ c^^Q-n m2 

Ji(e) , J2(g) 1 , 

C4, g g Ctf, g g 



Me) 



CiAl Q - n TO3 



Ccf, 



Q 



Uig) 



C1/.L Q - n 

c<j, g 



TO4 

Q ' 



(21) 



The spinor field X verifies the same equations apart from the fact that, due to its coupling to [see Eq. (^)], it is 
necessary to transform n ^ —n. 

As was originally found by Jackiw and Rossi (28j and Witten there are always n normalizable zero energy 
solutions of the Dirac operator in the vortex which allow fermions to propagate at the speed of light in the "— z" 
and say, directions, for the ^' and X fields, respectively. These solutions are found to be eigenvectors of the 

7°7'^ operator and are clearly obtained from the above equations by setting the consistency angular relationships 
TOi — 1 = 7714 + n and 7712 + 1 = TO3 + 77, those leading to the zero mode dispersion relation m = <^ cj = ±fc. Note 
that only one eigenstate of ^'^^'^ end up being normalizable for each kind of chiral coupling to the Higgs field, and 
thus the relevant dispersion relations reduce to uj = —k and w = fc, for the 5* and X zero modes, respectively [ p5[ . 

Such zero modes have a simple interpretation: since the Higgs field vanishes in the string core, the mass term Mjr 
in Eq. ( p^ vanishes too, and the fermions trapped in have zero mass. As a result, they propagate at the speed of 
light and they verify the dispersion relations ui = k on uj = —k. 



B. Massive trapped waves 



However, it is also possible a priori^ for the trapped fermions, to explore outer regions surrounding the string core 
where the Higgs field takes nonexactly vanishing values. In practice, this is achieved by means of a nonvanishing 
fermion angular momentum, which will lead to a nonvanishing effective mass to^ = uP' — k'^ ^ Q. For the ^ field, such 
massive solutions of the equations of motion ( po|) can only be obtained for four-dimensional solutions, in order to ease 
the zero mode constraint oj = ±k. The required angular consistency relations therefore read 

777 = 7771 = TO2 + 1 = TO3 + 77 = 7774 + 77 + 1. (22) 
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Similarly, the angular dependence of X field has to verify analogous conditions with the transformation n — > — n. It 
was previously shown numerically that the Abelian Higgs model with one Weyl spinor always admits such kind of 
normalizable solutions . In the following, massive solutions for Dirac spinors are numerically derived for our model 
and shown to exist for a wide range of fermion charges and coupling constants. 



1. Analytical considerations 



Some interesting analytical asymptotic behaviors of these modes have been previously studied |25 2^. In particular, 
there are only two degenerate normalizable eigensolutions of Eqs. ( po| ) at infinity. Since the Higgs field goes to its 
constant vacuum expectation value and the gauge coupling functions vanish, we found the eigensolutions to scale as 
exp (±rip), with 



n = 



m 



(23) 



First, note that in order to have decreasing solutions at infinity, the mass of the trapped modes m has to be less 
than the fermion vacuum mass rrif , as intuitively expected (for m > mf, one recovers the oscillating behaviour that 
is typical of free particle solutions). Moreover, from Cauchy theorem, two degrees of freedom can be set in order to 
keep only the two well defined solution at infinity. On the other hand, by looking at the power-law expansion of both 
system and solutions near the string core only two such solutions are also found to be normalizable. More 

precisely, normalizability of each eigensolution at g = leads to one condition on the winding numbers nii of each 
spinorial component ^i. Moreover, in order for the fermion field to be well defined by rotation around the string, 
each spinorial component with nonzero winding number nii has to vanish in the string core, and so behaves like a 
positive power of the radial distance to the core. The analytical expression of the eigensolutions near g = reads |2^] 



a3(al)^?""+l"l+2 

V a4(ai)f?-™+l"l+i J 

( aiQ"^ \ 
a2{ai)g^-^ 

a3{ai)g"'+\"\ 

V a4(al)^?™+l"l+^ / 




a2(ai)£<™+l"l-"+i 
«3(al)^'"-" 
Va4(al)^^"-"-^ / 

02(01)^""'+'"'+"+^ 



(24) 



The normalizability condition for the four eigensolutions can be summarized by 

sup (0, n) < m < inf (1, 1 + n), 



(25) 



and so, for any value of m there are only two conditions satisfied. However, from the consistency angular conditions 
on each spinorial components, only three pairs of solutions are acceptable near the string. Assuming n > 0, if m < 
then only the pair (si, S4) is both normalizable and well defined by rotation around the vortex, similarly for m > n + 1 
the relevant solutions are (32,33), whereas for 1 < m < n, they are (53,54). As a result, the two remaining degrees 
of freedom can be set to get only these pairs near the string core for a given value of m, but there is no reason that 
they should match with the two normalizable solutions at infinity. In order to realize this matching we have to fine 
tune another parameter which turns out to be the mass of the modes, m. As expected for bound states, this mass is 
therefore necessarily quantized. Note at this point that m = is, in such a procedure, nothing but a particular case 
of the general solution here presented. The three different pairs of well defined solutions at the origin suggest that 
there are three kinds of similar massive bound states in the vortex, according to the values of the winding number m. 
Intuitively, the more the field winds around the string, the farther the particle explores regions surrounding the core 
due to the higher values taken by its angular momentum, meaning the largest the extension of its wave function is, 
the more it acquires mass from coupling to a noncxactly vanishing Higgs field. As a result, the lowest massive modes 
will certainly be obtained from values of m which correspond to vanishing winding numbers m^. 
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2. Symmetries 



In the following, the equations of motion ( pO[ ) will be summarized in the form (5^)^ aj — 0, with implicit summation 
implied over repeated indices. 



The first symmetry is obtained from the complex conjugation of the equations of motion (20[). Since complex 
conjugation does not modify Eqs. (^0|), once the angular consistency relations ( p2| ) are set, there is an arbitrary 
complex phase in the choice of solutions, and it will be sufficient to look for real rescaled spinorial components 3^. 

There is another symmetry between the positive and negative energy solutions of the equations of motion ( pO| ) that 
may be useful. With the label e = ± for particle and antiparticle states, respectively, one has 

(5+)^ 5,^=0 ^ (5_)^ 5,„=0, (26) 

provided 

5.- =(7V)!5,,- (27) 

As a result, the negative energy solutions are obtained from the positive ones by the action of the '^'^^^ operator, thereby 
generalizing the properties of the zero modes which were precisely found as eigenstates of this operator |^,^,^ . 
The last symmetry concerns the gauge coupling functions fi. Under the transformations 

m —ffh — n + 1 ~ m, 

c^L ^ ^ ^^i'Rt (28) 

the gauge functions fi, in Eqs. (pO[), are simply swapped according to /i ^ and /2 ^ /a- As a result, for every 
5 solution found at given c^^ there is another solution a, with charge c^^ = — c^^ ^''^'^ winding number 

fh = n + 1 — m, namely 

(isig) = a2{g), a^ig) = ai{g). '-^^^ 

Note that the particular case c^^ = c^l = 2^/2 appears as a frontier separating two symmetric kinds of solutions 
with two differents winding numbers lying on both sides of m = (n + l)/2. As a result, the three different behaviors 
found above from normalization and angular consistency conditions seem to reduce to only two, since the domains 
where m < and m > n + 1 are actually connected by charge symmetry in relation to 

On the other hand, due to its coupling to the antivortex instead of the vortex, the equations of motion of the X 
field are simply obtained from Eqs. ( |2Q| ) by the transformations o^i — > c^^^^j^^ — > ^xl(r) ' — ^ ^i- The li are 

the winding numbers of the scaled X spinorial components, namely the (3i, and they verify the angular consistency 
relations (E2) with n replaced by —n as previously discussed. Let us introduce one more transformation on the ^' 
parameters, 

m —> fh — I + n, 

c^L ~^ '^ip-L — ^XRj (30) 

Naming Iji the scaled gauge coupling functions of the X spinor, the ^' ones are found to transform according to 
/i ~* ff3) /2 —^54, /a ~^ ffij and /4 — *■ 52- Thus, if the a are solutions of the ^' equations of motion ([2^), with m 
winding number and c^^ charge, then there exist (3 solutions for the X field with same mass m, provided I — m ^ n 
and = = ~ c^, and they read 

Pi{g) = a^ig), (32{g) -aiig), 

~ ^ ~ ^ (31) 

Psig) = ai{g), Piig) = -~a2{g). 

Owing to these symmetries, it is sufficient to study the ^ equations of motion (|^, for various values of the winding 
number m and for left-handed part charges, namely c^^j higher or equal than 
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3. Numerical methods 



In order to compute the relevant massive wave solutions for the 5* fermions on the string, it is necessary to solve 
first the vortex background. At zeroth order, neglecting the back reaction of the fermionic currents, and in terms of 
the dimensionless fields and parameters, the equations of motion for the string forming Higgs and gauge fields read, 
from Eq. 





IdH 


HQ 


de2 


Q dg 






1 dg 




dg^ 


g dg 





-H{H^ - 1), 



(32) 
(33) 



where mb = qcrpr] is the classical mass of the gauge boson. The solution of these equations is well known ||Tj 
and shown in Fig. ^ for a specific (assumed generic) set of parameters. 
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FIG. 1. The solutions of the field equations for the vortex background. The Higgs field, H, takes its vacuum expectation 
value at infinity and the gauge bosons condensate in the vortex. 



The system of Eqs. (gO|) being linear and involving only first order derivatives of the spinor components, a Runge- 
Kutta numerical method of integration has been used. However, as noted above, since we are interested only in 
normalizable solutions, it is more convenient to perform the resolution from an arbitrary cutoff at infinity, toward the 
string core. Let us introduce g^o, the cutoff value on the dimensionless radial distance. From the asymptotic form of 
Eqs. (po|) at infinity, and in order to suppress the exponential growth, the spinorial components cti have to verify 

~ f \ m ^ , m[ ^ 

Oii{goo) = ~t: a2(£'oo) + 7; a4(£'oo), (34) 

a3[Qoa) = -t: a2[goD) + -p. ai[goo)- (35) 

These conditions constrain two degrees of freedom, and another one is fixed by normalization of the wave functions 
at goo- As a result, only one free parameter can be used yet in order to achieve the matching between these well 
defined solutions and the two normalizable ones near the string core. It will be the case only for particular values of 
the mass m. Numerically, the matching is performed in two steps. First, by means of the last free parameter, one of 
the usually divergent component near the string core is made to vanish at g — 0. Obviously, this component is chosen 
among those having a nonzero winding number since, in order to be single valued by rotation around the vortex, it 
necessarily vanishes at the string core. Once it is performed, the last divergent component at = is regularized, 
its turn, by calculating the mass of the mode m leading to a convergent solution. For the range of model parameters 
previously defined, the numerical computations thus lead to the mass of the trapped wave solutions as well as their 
components as function of the radial distance to the string core ai{g). 
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4- Numerical results 



In what follows, the Higgs winding number is assumed fixed to the value n — 1, and the range of c^^ restricted to 
Cip-L ^ C0/2, the other case being derivable from the symmetric properties discussed above. 

The first results concern the "perturbative sector" where the fermion vacuum mass verifies mf < mh, or equivalently, 
for a smaller Yukawa coupling constant than the Higgs self-coupling, i.e., g < -s/A. In this case, for reasonable values 
of the fermion charges, i.e., of the same order of magnitude than the Higgs one c^^ ~ C0/2, only one normalizable 
massive bound state is found with null winding number m = 0. As a result, by means of transformations (p8|), there 
are also symmetric modes for c^,^ < with winding number m = 2. The dependency of the mode mass to with 

the fermion vacuum mass and charges (i.e., the coupling constants to Higgs and gauge fields) is plotted in Fig. ^ and 
Fig. 1^. The study has been also extended to the nonperturbative sector where this massive mode thus appears as the 
lowest massive bound state. First, it is found that the mass of the trapped mode always decreases with respect to 
the coupling constant, i.e, with the fermion vacuum mass TOf . Moreover, for small values of TOf/TOh, the derivative of 
the curve TO(TOf/TOh) vanishes near the origin (see Fig. |^). As a result, the mass modes in the full perturbative sector 
does not depend on the coupling constant to the Higgs field, at first order. On the other hand. Fig. ^ shows that the 
mass of the bound state hardy depends at all on its coupling with the gauge field (i.e., on the charges c^l) i'^ the 
nonperturbative sector, where all the curves have the same asymptotic behavior. Near the origin, the closest c^^ to 
the higher mode mass to is. In the particular limiting case c^-^^ ~ C0/2, there is no normalizable massive bound 
state, and as can be seen in Fig. ^, already for c^l/c^ = 2, the mode mass is close to TOf. It is not surprising since, as 
it was above noted, c^^ = C0/2 is a frontier between two kinds of solutions with different winding numbers, and thus, 
at this point, the "normalizable" winding numbers are not well defined. Note that this is only true if to 7^ (n + l)/2 
as it is the case here in the perturbative sector with n = 1 and to = 0. 
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FIG. 2. The mass of the lowest massive bound state, relative to the fermion vacuum mass, plotted as function of the coupling 
constant to the Higgs field, i.e., the fermion vacuum mass relative to the mass of the Higgs boson. 
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FIG. 3. The mass of the lowest massive bound state, relative to the fermion vacuum mass, plotted as function of the coupling 
constant to the Higgs field, for several values of the fermionic charges. The closest c^^^ is to 0^/2, the higher mode mass m is. 
In the extreme case c^^ ~ 2^/2, m ~ mf there is no longer normalizable massive bound state in the perturbative sector. 



The normalized scaled spinorial components a{g) have been plotted in Fig. ^ for the lowest massive bound state, 
with the normalization 

J gdgalai = 1. (36) 

The corresponding transverse probability density has also been plotted in Fig. ^. Note that the massive mode wave 
function is larger around the string rather on it, as expected for a nonvanishing angular momentum. 

0.15 T 1 I 1 




10 20 30 40 10 20 30 

P P 

FIG. 4. The transverse spinorial components of the ^I" field, as functions of the distance to the string core, for the lowest 
massive bound state. The transverse normalized probability density is also plotted and takes its maximum value nearby the 
string core, as expected for massive bound states exploring the neighborhood of the core by means of nonvanishing angular 
momentum. Note that for m = 0, one spinorial component behaves like a zero mode one, i.e., it condenses in the string core 
contrary to the others. 



The nonperturbative cases with mt > Wh, involve much more massive bound states. First, another mode appears 
in addition to the previous one, with the same winding number. Because of the fact that m decreases with mf [see 
Fig. ^, for higher values of mf , another mode comes into the normalizable mass range. Since normalizability at infinity 
requires m < mf , the number of massive modes increases with the value of mf . Moreover, there are also solutions 
involving all the other possible winding numbers. The evolution of the mass spectrum, for winding number m = 0, 
and with respect to the coupling constant to the Higgs and gauge fields is plotted in Fig. ^. The behavior of each 
mass is the same as that of the lowest mode previously studied, the new properties resulting only in the appearance of 
new states for higher values of the fermion vacuum mass mf , as found for two-dimensional Weyl spinors in Ref . pq| . 
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FIG. 5. The evolution of the mass spectrum for m — winding number, as function of the couphng constants. Each main 
branch represents one massive mode whereas the substructures show its evolution with respect to the fermion charge c^^^ . Five 
values of the fermion charge have been plotted, from c^^^ = 2 to c^^^ = 10, and the spectrum has been computed only in the 
nonperturbative sector, since only the lowest mode exists for lower values of mf/rrih (see Fig. ^). As expected, all the modes 
have mass fn decreasing with their coupling constant to the Higgs field. Moreover, the substructures show that, for sufficiently 
large values of mf/mh, the mode mass is a decreasing function of the charge c^^. However, note that this behavior can be 
inverted for some modes close to their appearance region, as it is the case for the second one. 

Physically, the additional massive modes at a given winding number can be interpreted as normalizable eigenstates 
of the angular momentum operator in the vortex potential, with higher eigenvalues. From Fig. |^ and Fig. ^, one 
can see that for each value of m, the lowest massive state is confined around the string with a transverse probability 
density showing only one peak whereas the higher massive modes have transverse probability density profiles with 
an increasing number of maxima, as can be seen in Fig. |^. In fact, as for the structure of atomic spectra, the two 
spatial degrees of freedom of the attractive potential certainly lead to two quantum numbers labeling the observable 
eigenstates, one of them being clearly to, and the other appearing through the number of zeros of the spinorial 
components, or, equivalently, the number of maxima of the associated transverse probability density. 
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FIG. 6. The transverse spinorial components of the field, plotted as functions of the distance to the string core, for the 
m = — 1 lowest massive bound state. The transverse normalized probability density is also plotted and vanishes in the string 
core. In this case, all components of the spinor wind around the string and the corresponding mode is thus centrifugally 
confined in a shell nearby the core, as expected for a nonzero angular momentum eigenstate. 
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FIG. 7. The components of the 'i' field as function of the radial distance, and the corresponding transverse probability 
densities. The curves have been computed for m = — 1 winding number, and for two additional incoming modes in the 
nonperturbative sector. As expected, interferences fringes appear from the nonzero angular momentum eigenvalues of these 
modes. 



The massive modes with higher winding numbers behave in the same way. However, they exist only for nonzero 
values of the coupling constant mf/mh, this one increasing with the value of the winding number m. The scaled 
spinorial components and the transverse normalized probability density of the lowest massive bound state with next 
TO = —1 winding number are plotted in Fig. |^. They are found to be normalizable for coupling constant mf/rrih ^ 0.5 
when c^^/cij, — 2, as can be seen in Fig. ^. Contrary to the m — lowest massive state, all spinorial components wind 
around the string, and the transverse probability of finding such a mode vanishes in the string core, as expected for 
a nonzero angular momentum eigenstate. Obviously, this is also true for all higher values of to, as for the to — —2 
massive mode which appears to be normalizable for mf/mi^ }^ 1.2. 




0.4 



mf/mh 

FIG. 8. The mass of the m = — 1 lowest massive bound state, relative to the fermion vacuum mass, with respect to the 
coupling constant to the Higgs field. Note that the mode does not exist in the perturbative sector. 



It is clear from the numerical results that the fermions can be trapped in the string in the form of massive bound 
states, for a wide range of model parameters. The only exception takes place for fermion charges close to the particular 
value = C0/2 where there is no normalizable massive bound state in the perturbative sector. Note, again, that all 
the previous results are also relevant for the massive modes with symmetric winding numbers m = n + l — m, as well 
as for the X spinor field and for the antiparticle states of both ^ and X fermions. 



IV. FOCK SPACE FOR MASSIVE MODES 

The existence of massive trapped waves requires that the quantum state space be enlarged to include them. 
For each normalizable mode with mass to, a two-dimensional Fock space can be constructed by spinor field expansion 
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over the relevant massive plane waves. The full quantum theory can therefore be obtained from tensorial product of 
the different Fock spaces belonging to their own mass representation, together with the Fock space associated with 
the zero modes As a first step toward a full theory, we will only consider the plane waves associated with one 
massive mode of mass m. 



A. Quantum field operators 



Quantization is performed through the canonical procedure by defining creation and annihilation operators satisfying 
anticommutation rules. However, the particular structure of the trapped massive waves yields relationships between 
longitudinal quantum operators with nontrivial transverse dependencies. 



1. Fourier transform 



In the previous section, it was shown that the fermions could propagate along the string direction with given mass 
m belonging to the spectrum. From Eq. ( p^ ) , setting 



with 



(37) 



(38) 



and using the symmetry properties shown in Sec. (IIIB2), the transverse parts of the massive trapped waves for 
particle and antiparticle states can be written as 



u^{k,x±_) 



with the notations 



/ Vt^r+i;ai(r)e-™^ 
^^^^7^a2(r)e-^('"-l)^ 
V^w^a3(r) e-»(™-")« 
V iVwna4(r)e-*(™-"-i)^ / 



\ 



V iV'ir+fca4(r)e-*('"-"-i)'' / 



a;_L = (?■, ^), and a — — 



'2tt 



(39) 



(40) 



Contrary to the zero mode case, fermions can now propagate in both directions of the string, so that the momentum 
k of the massive waves can take positive and negative values. As a result, the ^' field can be Fourier expanded over 
positive and negative energy states as 



dk 
2tt2uj 



b\k) u{k, xi) e*M-'=-) + d{k) v{k, XA_) e-^i-'i-i'^) 



(41) 



where the subscripts have been omitted. The normalization convention of the Fourier transform is chosen as in the 
zero modes case [g5[, i.e., 

j dze*('=''=')^ = 27r5(fc - k'). (42) 
Obviously, the X field verifies similar relations with the transformation n —f —n, as was noted previously. 



2. Commutation relations 



In order to express the Fourier coefficients b{k) and d^(fc) as function of the spinor field let us introduce another 
unit spinors 



u{k,x±) = 



\ 



I -v/wTfcS3(?-)e~™® 
iVtr^a4(r)e-*(™-i)'' 
^/w^al(r) e-^('"-")» 
V «Var+fca2(r) e-H'"-«-i)e j 



v{k,x±) = 



( + ka-i{r) e" 



imG 



(43) 
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They clearly verify u = ^^j^v and from Eq, 

u'^{k)u(k) = vHk)v{k) = 2ujV{r) 
u{k)^v{-k) = d^k)u{-k) = 0, 



U,(-lc\ = ^Tt(''^^»■^ - r, (44) 



where the dependency with respect to transverse coordinates have been omitted in order to simphfy the notation, 
and where we introduced the function 

77(r) = ai {r)a3 (r) + a2 (r)a4 (r) . (45) 
From Eqs. (^), (^), and (0), the Fourier coefficients are found to be functions of the '5 field, and read 



b''{k) ^ Jf j rdrdedze-'^'^^-'^''^u\k,xi_)-^, 
d{k) = J rdrd0dze'^'^'-'''^v\k,x^)^, 



(46) 



where we have defined the normalization factor 

A/'= / rdrdev{r)= I gdgli{g) (47) 



Similarly, the expansion of the ^E"^ field on the same positive and negative energy solutions leads to the definition of 
its Fourier coefficients, namely b{k) and S{k). From Eqs. ( ^ ) and (^), they can also be expressed as functions of 
^\ and verify 

b{k) [b\k)]\ and d^k) = [d{k)]K (48) 

In order to perform a canonical quantization along the string world sheet, let us postulate the anticommutation 
rules, at equal times, between the spinor fields 

{vl/,(i,x),vl/«(t,f')} ^6{z^z'){ry^{x^,x'^), (49) 
where F" is a matrix with respect to spinor components whose utility will be justified later, and which reads 

r"(xi,a;l)-^(c^I-fc7V) [u{k,x±)u\k,x'j^) + v{k,x^)v\k,x'^)] . (50) 

Note that F'^ does not depend on uj and k. Explicit calculations show that the first terms involving lo and k are mixed 
with u{k) and v{k), and yield Lorentz invariant quantities, such as m. Moreover, F" has the following orthonormal- 
ization properties 



u^{k,xi_)T°{xj_,x'^)u{k,x'^) = 2ujV{r)v{r') 
u^{k,x±^)T'^{x_i_,x'^)v{~k,x'j_) — 0, 



(51) 



and similar relationships are obtained for v by swapping u and v. 

The anticommutation rules for the Fourier coefficients are immediately obtained from Eqs. (^), (^), (^), and 
Eq. (^), using the properties of F° in Eq. (|5ll), and read 

{b{k),b^k')} = {d(fc),dt(/c')} = 2n2ujS{k~ k'), (52) 

with all the other anticommutators vanishing. As a result, the Fourier coefficients b^ and d^ behave as well defined cre- 
ation operators, whereas their complex conjugates, b and d, act as annihilation operators of a particle and antiparticle 
massive state, respectively. 

In order to verify the microcausality of the theory and to justify, a posteriori, the ansatz of Eq. (^9|), let us derive 
the anticommutator between the quantum field operators 5" and at any time. The ^E" expansion in Eq. (|l]) and 
its complex conjugate yield 
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{vl/,(x),M't.(x')} 



dfcdfc' 



{&t(fc),6(fc')}u»(fc,a;i)utJ'(fc',a;'L)e*[('^*-'^'*''-(*=^-'='^')] + 



(27r)22w2w' 

{d(fc),dt(fc')}„,.(fc,a;i)wtJ'(fc',a;l)c-'[('^*-'^'*')-(*=^-'='^')l 



(53) 



Using Eq. (p2|), this equation simplifies to involve tensorial products of unit spinors evaluated at the same momentum. 
It is therefore convenient to introduce two additional matrices, namely T^{x±,x'j_) and A4{x±,x'j_), which verify 



u{k, x±)u^ {k, x'j_) — ujT'~'{x±,x'^) — kr^{xj_,x'^) — Ai{xj_,x'^) 
v{k,x±)v''{k,x'^) = ujT°{x±,x'j_) - kr^{x±,x'^) + M{x±,x'^). 

From Eq. (|39|), these matrices are simply related to r° by 

r^{x^,x'^) = rO(:E^,a;l)7V, 
Mix^,x'J r'>ix^,x'JMdix'^h°, 

where Md{x±) is the diagonal matrix 



(54) 



^3 (^) 

Mdix±) — mDiag ( _ , , e 



aijr) ^^,„e aijr) ^,„g 0^2 (0 ^,„e 
^ai(r)'' 'a2(r) ' a3(r) ' a4(r) 

From Eqs. ( p^ ) and (^5[), the anticommutator ( p3| ) reduces to 

{^P(x),vl.t(x')} ^ [rO(:E^,a;'^)*ao+r3(a;^,xl)»53+A^(xi,xi)] zA(a;|| - rcj,), 

where a;|| = (i, z), and A is the well-known Pauli Jordan function which reads 



iA(a;|| — ) 



dfc 



(55) 



(56) 



(57) 



(58) 



and vanishes outside the light cone. As a result, the quantum fields indeed respect microcausality along the string. 
The matrices F^ appear as the analogues of the matrices 7'' for the Dirac spinors living in the vortex. The two- 
dimensional quantization along the string is thus not independent of the transverse structure. It is all the more so 
manifest in the anticommutator expression between ^' and from Eq. (57), and using Eq. (pa), one gets 



{*(x),^'(x')} ^T°ix^,x'^){zj°do+ij'd3+Md{x^)) iA(a;|| -xj|). 



(59) 



The matrix T° now appears clearly as a local transverse normalization of the longitudinal quantum field operators. 
Note that the mass term also depends on the transverse coordinates due to the nontrivial profile of the Higgs field 
around the string. Moreover, setting t = t' in Eq. (^, leads to the postulated anticommutator at equal times (|4S|), 
and therefore justifies the introduction of the F" term. 



3. Fock states 



In the following, \V) will design a Fock state constructed by applying creation operators associated with a massive 
mode m, on the relevant string vacuum. Such a state was similarly defined for zero modes in Ref. [^ . From the 
anticommutators (p^), a massive \E' state with momentum k is now normalized according to 

{k'\k) ^2Tr2ujS{k-k'). (60) 

Similarly, it will turn out to be convenient to derive the average of the occupation number operator since it will appear 
in the derivation of the equation of state. From Eq. (p2|), and for a ^ massive mode, it reads 



where the summation runs over all 5* massive particle states present in the relevant m Fock state {V), and L is the 
physical string length, coming from the 5(0) regularization by means of Eq. (42). 
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B. Stress tensor and Hamiltonian 



The classical stress tensor can immediately be derived from variation of the full Lagrangian (|l|) with respect to the 
metric, and the ^' fermionic part thus reads pq| 



(62) 



1. Hamiltonian 



The quantum operators associated with the classically conserved charges can be obtained by replacing the classical 
fields by their quantum forms involving creation and annihilation operators. In this way, the Hamiltonian appears, 
from Noether theorem, as the charge associated with the time component of the energy momentum tensor 

T^* = i^fdt"^ - i (a**) 7"*. (63) 



Using Eqs. ( |39| ) and (41) in the previous equation, and performing a spatial integration, the Hamiltonian operator 
reads, after some algebra, 

/Ah r 
^ / d?x^ [-h{k)h\k) + d\k)d{k)\ uik, x^)-f%ik, x^). (64) 

In order to simplify this expression, let us introduce the parameters 

Ex = ^2 + ^3, and Sy — af + a\. (65) 
From Eq. (p^), the Hamiltonian now reads 

/Ah 
^ [-b{k)h\k) + d)[k)d{k)] [{u^ - k)\\Y.l\\ + (c + k)\\Y.\\\] , (66) 

with 

' rdrdOT: ^ j QdQY^^g). (67) 

Analogous relations also hold for the X field. It is interesting to note that Eq. (|6^ ) generalizes the expression previously 
derived in the zero modes case [ p5| , being found again by setting lo = —k for the ^' zero modes, 01 to = k for the X 
ones. The normal ordered Hamiltonian is obtained if one uses the anticommuting normal ordered product for creation 
and annihilation operators, i.e., 

■■K--- j ^[b\mk)+d\k)d{k)] [{u:-k)\\Y.l\\ + {Lo + k)\\Y.\\\]. (68) 

Since uj > k, this Hamiltonian is always positive definite and is thus well defined. Note that, as in the zero mode case, 
such a prescription overlooks the energy density difference between the vacuum on the string and the usual one, but 
it can be shown to vary as 1/L^ and therefore goes to zero in the infinite string limit p^,^Q . 



2. Effective stress tensor 

In order to make contact with the macroscopic formalism p^ , it is necessary to express the classically observable 
quantities with no explicit dependence in the microscopic structure. The relevant two-dimensional fermionic energy 
momentum tensor can be identified with the full one in Eq. ( |6^ ) , once the transverse coordinates have been integrated 
over. Due to the cylindrical symmetry around the string direction, z say, all nondiagonal components, in a Cartesian 
basis, involving a transverse coordinate vanish after integration. Moreover, since the fermion fields are normalizable 
in the transverse plane, the diagonal terms, and T^^, are also well defined, and by means of local transverse 
stress tensor conservation, the integrated diagonal components, T^^ and T^^, also vanish As expected, the only 

relevant terms in the macroscopic formalism are thus T^^ with a, f3 E {t, z}, i.e., the ones that live only in the string 
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world sheet. On the other hand, the macroscopic Umit of the involved quantum operators is simply obtained by taking 
their average over the relevant Fock states. 

Replacing the quantum fields in Eq. (^) by their expansion (^H), and using Eq. (^), one gets the quantum 
expression of the energy momentum tensor. Averaging the relevant components T^^ in the Fock state jP), by means 
of Eqs. (11) and (|l]), one obtains 



1 

L 

1 
L 



N, 



—u{ki)-i^u{ki) + —uikjh^u{kj) 



u{ki)^ u{ki) H u(ki)^ u(ki) 



+ {i,N^,)^{3,N^) , 



(69) 
(70) 
(71) 



where the i summations run over the particle states with momentum ki involved in the Fock state jP), while the 
j summations take care of the antiparticle states with momentum kj , all with mass m. In order to simplify the 
notation, the transverse dependence of the unit spinors have not been written, and the averaged operators stand for 



(72) 



Similarly, the same relationships can be derived for the X field, by replacing the ^ unit spinors by the X ones with 
the correct angular dependence, and certainly, in another mass representation m^. Once the transverse coordinates 
have been integrated over, Eqs. (|69|), ([70|), and ([7l|), lead to the two-dimensional 4* stress tensor 



/ 



E^,^\\^'y\\ + E^^\\i:i\ 



E. 



P 



ip-p 



E. 



'Ip-p Pip-p 



p^p\\j:1\\ - p,i,j\j:1 



(73) 



with the notations 



E^p — 



1 

L 



N,, 



V — (wi + fc,) + V — {iOj + kj] 



E. 



Ip-p 



1 
L 



P-4'V ^ T 



k ""^ k 



(74) 



Recall that the full effective energy momentum tensor also involves the Higgs and gauge fields of the vortex background. 
Since they essentially describe a Goto-Nambu string [R2| , their transverse integration yields a traceless diagonal tensor 



J rdrdO (T** +Tl')=- j r dr de* (r|^ + TjJ") = . 



(75) 



Note that the full stress tensor may also involve several massive ^I' and X states, with different masses belonging to 
the spectrum. In this case, there will be as many additional terms in the form of Eq. (73), as different massive states 
there are in the chosen Fock states. 



C. Fermionic currents 



The quantum current operators can be derived from their classical expressions (^4|) by using Eq. (|l|), while the 
corresponding conserved charges are obtained from their spatial integration. By means of Eq. (^), the current 
operator, averaged in the relevant Fock state \V), reads 
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' 2 L 
' 2 L 



I,): ^-^ I,): 



I,): ^—^ I,): 



(76) 



with a e {t, 2:}, and once again, the sums run over or X, particle and antiparticle states. The Ui are the unit 
spinors associated with the field dealt with. Concerning the transverse components, due to the properties of the unit 
spinors u and v in Eq. P9), only the orthoradial one does not vanish and reads 



]jr ■.)v = -q 



q- 



2 L 

c^R - cjfl 1 



— g N:r — g 



Nj^ _ Q N:f — g 



(77) 



whereas (: :) = due to the bound state nature of the studied currents. As expected, the gauge charges carried by 
each trapped fermion in the form of massive mode, generate only macroscopic charge and current densities along the 
string, as was the case for the zero modes However, the nonvanishing orthoradial component shows that the local 
charges also wind around the string while propagating in the z direction, as suggested by the above numerical studies. 
However, this component will be no longer relevant in the macroscopic formalism, since it vanishes in a Cartesian 
basis, once the transverse coordinates have been integrated over. 

Nevertheless, this nonzero angular momentum of the massive modes is found to generate new properties for the 
longitudinal currents. Let us focus on the vectorial gauge currents generated by one exitation state, with energy to 
and momentum fc, of a ^' massive mode, m say. From Eq. (76), using Eqs. ( ^9|) and (|65|), the world sheet vectorial 
charge current reads 



3%v '■) 



' 2 L 

c^R + £ 
' 2 L 



1 



k \ —2 



1 - 



k \ —2 



UJ 



-'X 



UJ 



(78) 
(79) 



where e = ±1 stands a one particle or antiparticle exitation state. Now, even setting fc = in the previous equations 
yields a nonvanishing spatial current. Physically, it can be simply interpreted as an anomalous magneticlike moment of 
the considered massive mode in its rest frame. Examining Eq. ( [t^ shows that it could be null only if 'Syi^) = ^x('')j 
which is generally not satisfied due to the particular shapes of massive spinor components trapped in the string (see 
Sec. [H). These ones being associated with nonzero winding numbers, it is therefore not surprising that, even for a 



massive stationary state along the string, the nonvanishing charge angular momentum around the string generates such 
additional magneticlike moment. Note that it does not concern the zero modes, first because they precisely involve 
vanishing winding numbers ||25[| , and then because for them, there is no defined rest frame due to their vanishing 
mass. Obviously, this property can be generalized for the axial part of the current, and thus is also valid for the total 
current of any massive spinor field trapped in the string. 

All the above construction of the Fock space and the derivation of the quantum operators associated with the 
energy momentum tensor and gauge currents remains valid for each 5* and X massive mode. More precisely, the other 
masses belonging to the ^ spectrum verify analogous relationships provided m is replaced by the relevant one, as for 
the unit spinors. In addition, the X massive states require to transform n — > — n in Eq. (56), due to their coupling to 
the antivortex. At this stage, the averaged values of the stress tensor and currents have been obtained, and therefore 
allow the derivation of an equation of state, once the Fock states are known. 



V. EQUATION OF STATE 

The energy per unit length and tension in a given Fock state \V) are basically the eigenvalues associated with 
timelike and spacelike eigenvectors of the effective two-dimensional full stress tensor. Obviously this one includes the 
classical Goto-Nambu term resulting of the string forming Higgs and gauge fields [see Eq. (|75|)], with the fermionic 



17 



part generated by the massive currents [see Eq. (|7^)]. Moreover, in order to describe the string by an adequate 
macroscopic formahsm, it is necessary to choose a quantum statistics for the relevant Fock states, and for energy 
scales far below the ones where the string was formed, it is reasonable to consider a Fermi-Dirac distribution at zero 
temperature |^,^. In the following, the equation of state is first derived for the lowest massive modes associated 
with the ^' and X field, and simplified in the zcro-temperaturc and infinite string limit. As a second step, these 
derivations are generalized to any number and kind of trapped fcrmionic mode. 



A. Lowest massive modes 



1. Energy per unit length and tension 



In this section, we will only consider the lowest massive modes belonging to the and X mass spectrums, with 
masses and rn^, respectively. From Eqs. ( |73| ) and (75), the full effective energy momentum tensor reads 



(80) 



where takes the same form as t'^ in Eqs. @ and (0) once the 5" relevant parameters have been replaced 
by the X ones. In the preferred frame where the stress tensor is diagonal, we can identify its timelike and spacelike 
eigenvalues with energy per unit length U and tension T. Upon using Eqs. (73), (|75|), and (80), these read 



Uv = Af 2 + 



E:f-P:f, 



J2 {E:f + P:f)\\^ 



1/2 



(81) 



for the energy per unit length, and 



Tv = Ar 



Ejr + Pr, 



1 1/2 



E {E:,-Pr)\\^l 



(82) 



for the tension. It is interesting to note first that U-p+T-p ^ 2M^, and thus the fixed trace equation of state previously 
found for zero modes is no longer verified by massive modes, as expected since they are no longer eigenstates 

of the 7^7'^ operator. Moreover, the expression of energy density and tension does not seem to involve the conserved 
charge current magnitude, which played the role of a state parameter in the case of a scalar condensate in a cosmic 
string ||l^,|l5|. In fact, as it was the case at zeroth order for the zero modes the charge currents are only involved 
in the stress tensor through their coupling to the gauge field [see Eq. (^2|)]. At zeroth order, when the back reaction 
is neglected, the only nonvanishing component of the gauge field is Bg, and it therefore couples only with j^, which 
vanishes once the transverse coordinates have been integrated over. As a result, it is not surprising that the equation 
of state does not involve the fermionic currents without back reaction. As a result, it is more natural from quantization 
to define the occupation numbers of the involved species as state parameters. 



2. Zero -temperature and infinite string limit 

Assuming a Fermi-Dirac distribution at zero temperature for the exitation states, the sums involved in Eq. ( |7^ ) 
run over the successive values of the allowed momentum ki until the Fermi level of the considered species is reached. 
With periodic boundary conditions on spinor fields, the allowed momentum exitation values are discretized according 
to 

27r 

kn = -jjn, (83) 
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where n is an integer, playing the role of a quantum exitation number. As a result, in the relevant m representation 
of each field, the exitation energies uii are also discrete according to Eq. (|3^), and for the field, the parameters 
and in Eq. (|7j) simplify to sums of radical function of n, with n running from the vacuum to the last filled state. 
In order to express them as explicit functions of the relevant Fermi level, it is convenient to consider the infinite string 
limit L — » DO. In this limit one gets 



n + 

1 , 1 f ^^i> 



dkfik), 



(84) 



where — N^^ / L are the ^ up and down mover densities, , standing for the number of ^ particle moving in 
the +z or —z directions, respectively. Note that the total number of particles of this kind is thus iV^ — + + 1 
since there is the additional rest state obtained for k ~ 0. After some algebra in Eq. (^^, using Eq. (^^, the 
parameters, for the 4* field, read 



E,h — - — 



4^ 



^ An 



P^ - P4, 



where stands for the dimensionless ^' mover density 



27r 



In 



- In 



Pi> = —P^, 
m 



+ Pi + P^ 



+ Pi, +Pi, 



+ Pi + P^ 



+ P4> +Pi, 



(85) 



(86) 



(87) 



while is defined in the same way for the antiparticle states. Similarly, the two other parameters and 
read 



E,i 



-Pi, +A0 + 



m 
An 

( 



pU^ + pT +Pi.^^ + Pt +1^^ {j^ + pf+pl) (^i + pZ+Pv.) ) 



Pi^y^ + pi +Pv>V + 



In 



^^ + pX+pI) (\/i + p7 + Pv.) ) 



^ P^ 



(89) 



Note that these parameters depend differently on the up and down mover densities as expected for chiral coupling 
of the fermions to the string forming Higgs field. Recall that in the massless case the zero modes associated with 
the ^ and X fields can only propagate in the —z and +z direction respectively |^,^||^p8|. The same relationships 

also hold for the X field by using the relevant dimensionless mover densities and p^ . Although the equation of 
state can be derived as a function of these four parameters for each fermion field T , it is convenient at this stage to 
perform some physical simplifications. Contrary to the zero mode case, the coupling between massive particles and 
antiparticles of the same species T does not vanish along the string. As a result, it is reasonable to consider that the 
only kind surviving at zero temperature corresponds to the one which was in excess in the plasma in which the string 
was formed during the phase transition. On the other hand, the energetically favored distribution at zero temperature 
involves necessarily the same number of T "up" and "down" movers, each filling the accessible states living on each 
branch of the mass hyperbola (see Fig. ^). As a result, in the considered energy scale, it seems reasonable to consider 
only one state parameter per mass instead of the four initially introduced by quantization, namely pjr = 'pjr = p^, 
for a plasma dominated by particles, say. 
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FIG. 9. The filling of massive states trapped in the string, as expected in the zero-temperature limit, for a particular mass 
and for one species, ^' or A" say. All antiparticles have disappeared by annihilation with particles during cooling, and the 
interactions between particles moving in opposite directions, as their coupling to the gauge field, lead to the energetically 
favored configuration with same number of up and down movers. Obviously, the Fermi level is necessarily below the vacuum 
mass of the relevant fermion. 



Setting these simplifications in Eqs. ( p5D to (|89|), by means of Eqs. ( plD 
associated with the lowest massive modes now read 



and the energy density and the tension 
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The sum runs over the two lowest massive bound states, each one being associated to the two fermion fields trapped 
in the vortex, namely VE" and X, and have and masses, respectively. As a result, the equation of state involves 
two independent parameters, and p-^, in the zero-temperature and infinite string limit. The energy per unit length 
and the tension have been plotted in Fig. for the lowest massive modes in the nonpcrturbative sector. The curves 
are essentially the same in the perturbative sector, but the variations around the Goto-Nambu case [/ = T = A/^ are 
much more damped. For reasonable values of the transverse normalizations, e.g., ||S^|| ~ HSyll ~ 0.5, and for small 
values of the dimensionless parameters p^ and p^^, the energy density is found to grow linearly with p^ and p^^, whereas 
the tension varies quadratically. As can be seen in Eq. (Bffl, due to the minus sign in T, all linear terms in p vanish 
near origin, whereas it is not the case for the energy density. However, for higher values of the densities, the quadratic 
terms dominate and both energy density and tension end up being quadratic functions of p. On the other hand, 
according to the macroscopic formalism [p^ , the string becomes unstable with respect to transverse perturbations 
when the tension takes on negative values, as in Fig. |l^ for high densities. Moreover, the decrease of the tension is 
more damped in the perturbative sector, and the negative values cannot actually be reached for acceptable values 
of p, i.e., p < TOf/m. As a result, the rapid decrease of the tension with respect to the fermion densities constrains 
the nonpcrturbative sector where the string is able to carry massive fermionic currents. For each mass, the higher 
acceptable value of the p ensuring transverse normalizability is roughly mf/m, and from Eq. (|90|), the tension becomes 
negative at this density for mj ^ 47rM^. Much higher values of Wf will thus yield to empty massive states. 
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of the dimensionless effective densities of the two fermion fields, and p^- The parameters have been chosen in the nonpertur- 
bative sector, with m{/mii ~ 3 and fn^ ~ ~ 0.6mf. Note the linear variation of the energy density near the origin whereas 
the tension varies quadratically. Moreover, in the allowed range for fermion densities, i.e., less than the fermion vacuum mass, 
the tension vanishes and the string becomes unstable with respect to transverse perturbations. 



As previously noted, the energy density and tension for massive modes no longer verify the fixed trace equation of 
state found with the zero modes alone. As a result, the longitudinal perturbation propagation speed = —dT/dU 
is no longer equal to the speed of light, and it is even no longer well defined since the equation of state involves more 
than one state parameter. A necessary condition for longitudinal stability can nevertheless be stated by verifying that 
all the perturbation propagation speeds —{dT / dp) / {dU / d])) obtained from variation of only one state parameter are 
positive and less than the speed of light. The longitudinal and transverse perturbation propagation speeds, Cl and 
— T/U, respectively, have been plotted in Fig. 11^ in the case where there is only one species trapped as massive 
mode, ^! OT X say. It is interesting to note that there is a transition between a supersonic regime obtained at low 
fermion density, and a subsonic at high fermion density. Moreover, the transition density between the two regimes 
is all the more so high as the coupling constant mf/rrih is weak. It is not surprising to recover such zero- mode-like 
subsonic behavior [ES 33 1 for densities much higher than the rest mass, since in these cases the ultrarelativistic limit 
applies. On the other hand, since the mass of the massive mode decreases with the coupling constant as in Fig. g the 
transition will occur earlier in the nonperturbative sector, as can be seen in Fig. |l^. Note that the subsonic region 
is also limited by the maximum allowed values of the massive fermion densities, i.e, ~ rrif/rn, and the regions of 
transverse instabilities where becomes negative. Inclusion of the other species does not change significantly these 
behaviors, the main effect being to lower with respect to the other fermion density, as can be seen in Fig. nO. 
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FIG. 11. The squared longitudinal and transverse perturbations propagation speeds for one massive species only, plotted as 
functions of the dimensionless fermion density p, for two values of the coupling constant mf/mh. Note the transition between 
subsonic and supersonic behaviors takes place at a cross density, px say, which decreases with the coupling constant (see 
Fig. |l|). 
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FIG. 12. The cross dimensionless density px , i.e., the dimensionless fermion density for which the transverse and longitudinal 
perturbation propagation speeds are equal, plotted as function of the coupling constant mf/mh, for one massive species only. 
The dashed curve shows the maximum allowed values mt jm of the massive fermion density ensuring transverse normalizability. 
The transition from supersonic regime to subsonic can thus occurs only in the nonperturbative sector, below this frontier. 



The current magnitude can also be derived from the averaged current operators in the zero temperature and infinite 
string hmit. From Eq. ([T^), using Eqs. (^) and (|6l|), once the transverse coordinates have been integrated over, the 
world sheet components read 



-3 ^ ^ (^2) 

(j > = -Xl~('?^^Rll"i""2ll+9c^Ll|a4-a3ll) 



The current magnitude = (j")^ — (j^)^ therefore reads 



= 4 1 y: ^^Y^ j (^E ^^x^ p^) ' (93) 



where Fy^jr and Fyt denote the transverse effective charges 

Fyt = qc:F„, 1 1 a? 1 1 + QCj^l 1 1 "4 1 1 
Fxj^ = 'JCjFr 11^211 + gc^L II "si I' 



(94) 



already introduced for the zero mode currents |2q| . In the case of one massive species only, the charge current 
magnitude simplifies to 

C2 = A^F^FYfP, (95) 

and thus the sign of is only given by the sign of FxFy, which is generally positive for reasonable values of the 
transverse normahzations, e.g., ||Sx|| ~ H^yll ^ 0.5. As a result, the charge current generated by only one massive 
species is always timelike [26| , contrary to the zero mode charge current which was found to be possibly timelike, but 
also spacelike |25|, owing to the allowed exitations of antiparticle zero mode states. As noted above, the antiparticle 
states cannot exist for massive modes due to the nonvanishing cross section along the string between massive particles 
and antiparticles. Moreover, and as it was the case for the zero modes, unless there is only one massive species 
trapped in the string, the magnitude of the charge current is not a sufficient state parameter, contrary to the bosonic 
current-carrier case uM. 



B. General case 



From the numerical approach in Sec. (IIIB4), as soon as the nonperturbative sectors are considered, additional 
massive bound states become relevant, and it is reasonable to consider that, in the zero-temperature limit, all these 
accessible massive states will be also be filled. Moreover, the complete description of the string state also requires the 
inclusion of the zero modes in addition to the massive ones. 
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1. Full stress tensor 



The effective two-dimensional energy momentum tensor, involving all trapped modes in the cosmic string, can be 
obtained from Eq. (|o|) by replacing the sum over the two lowest massive modes with the sum over all the accessible 
masses, plus the zero mode terms previously derived in Ref. [ p5| . In the preferred frame where the stress tensor is 
diagonal, after some algebra, the energy density and tension therefore read 

1 
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TT 
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1/2 



(96) 
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(97) 



The sums run over all accessible massive bound states t with masses ffij^i of each fermion i.e., 4" and X. The 
additional parameters TZ^ and TZ^ are the particle densities trapped in the string in the form of zero modes, for the 
X and VE" field, respectively, with same notation as in Ref. ||2^. Note that the zero mode contribution can also be 
obtained from the null mass limit in Eq. (po|). As a result, the full expression of energy per unit length and tension 
seems to involve as many state parameters as trapped modes in the string. 



2. Equation of state 



As for the lowest massive modes, it is convenient to perform some approximations owing to the energetically 
favored filling of the involved states, in the zero-temperature limit. In particular, it is reasonable to consider that the 
nonvanishing cross sections between massive modes, and between zero modes and massive modes, lead to the filling 
of all the accessible states with energy lower than a Fermi energy, £jr say, for each fermion field J-. As a result, the 
energetically favored filling takes place by successive jumps from the lower masses to the highest ones, until the last 
mass hyperbola with mjr^ ^ Ejr is reached. Obviously, this filling begins with the zero modes, next with the lowest 
massive modes and so on. On the other hand, only the particle s tates a re assumed to be relevant because of the 
assumed annihilation of the antiparticle states, as discussed in Sec. ( V A 2 ). As a result, the Fermi levels, vj: say, can 
be defined through the zero modes filling only, as the line densities of zero mode exitations trapped in the string (see 
Fig. nj), and thus play the role of state parameters. 
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FIG. 13. The accessible states, for the X fermions, in the zero-temperature limit. The zero modes are represented by the 
chiral line uj — k, while the massive modes appear as mass hyperbolae. The Fermi levels are therefore dependent of the 
considered energy scale £, since the filling is performed by successive jumps from the zero modes to the massive ones, with 
fne ~ f . As a result, under these approximations, each trapped species leads to only one state parameter which can be identified 
with the Fermi level of the zero mode exitations, namely u = £/2tt. Note that the antiparticle states have not been represented 
due to their assumed annihilation. 



According to the so-defined state parameters, the massive exitation densities pjf^, in Eqs. (96) and (P7|), reduce to 



2tt 
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(98) 



with function is tlie Heavyside step function, as expected for energy scales less than the rest mass of the considered 
massive mode. The zero mode density simply reads 

Tljr ^ Vjr^ (99) 

for zero mode particle states alone. From Eqs. (|9|) and (p9|), and the definition of the dimensionless densities in 
Eq. (|^), 



27r 



(100) 



the energy per unit length and the tension now depend expHcitly of the two state parameters only, namely v^, and 
. By means of Eq. (pq) , the energy density reads 
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while the tension is obtained from Eq. (|97|), 
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The full energy per unit length and tension have been plotted in Fig. |lj for a configuration including two massive 
bound states, in addition to the zero mode ones. Due to the zero-temperature limit, for d ensiti es s maller than the 
first accessible mass, the Heavyside functions in Eq. ( p8|) vanish, as a result, from Eqs. (101) and (102) the fixed trace 
equation of state is recovered |23| with 



T 



(103) 



Once the first mass hyperbola is reached, the behaviors of the energy per unit length and tension are clearly modified 
and become very rapidly dominated by the mass terms, and, as found for the lowest massive modes alone, the energy 
density begins to grow linearly with respect to the state parameters, whereas the tension decreases quadratically. 
Actually, the plotted curves in Fig. |lj show slope discontinuities each time the phase space is enlarged due to the 
income of accessible massive bound states (see Fig. O). 
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FIG. 14. The energy per unit length and the tension plotted as function of the two state parameters, i.e., the zero mode 

densities, in the nonperturbative sector rrif/mh ~ 4. Two additional massive bound states have been considered with respective 

masses rn/mt ~ 0.4 and rn/rrif ~ 0.6. In the zero-temperature limit, the filling of the accessible states is performed by 

successive jumps as soon as the Fermi level reaches one mass hyperbola (see Fig. O). As a result, for the lowest values of the 

state parameters, only the zero modes are relevant and the fixed trace equation of state, U + T = 2M^, is verified, then the 

first and second massive modes are successively reached and become rapidly dominant. As can be seen near the origin, the 

smooth variations induced by the zero modes appear completely negligible compared to the massive ones. In the perturbative 

sectors, these behaviors are essentially the same, but the induced variations of the density energy and the tension are all the 

more small. 



On the other hand, it is reasonable to expect competition between the subsonic regimes induced by zero mode 
currents, or ultrarelativistic massive modes, and the supersonic ones coming from massive currents. In all cases, when 
the state parameters remain small, only the chiral massless states are accessible and the regime is obviously subsonic, 
as can be seen in Fig. ^ However, the massive mode filling modifies radicaUy this behavior, and as found for the 
lowest massive modes alone, as soon as a mass hyperbola is reached, the longitudinal perturbations propagation speed 
falls drastically and ends up being less than the transverse perturbation velocity. There is a rapid transition from 
the subsonic to the supersonic regime. For higher densities i^, the behavior depends on the coupling constant. More 
precisely, in the nonperturbative sector, the ultrarelativistic limit can be applied before the energy scales reach the 
fermion vacuum masses, and thus the subsonic regime is recovered, whereas it is not the case in the perturbative 
sector, as can be seen in Fig. 11^. 
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FIG. 15. The squared longitudinal and transverse perturbation propagation speeds for a spectrum involving two massive 
states in addition to the zero mode, plotted as functions of the state parameter ujr of one species, the other being fixed to 
a particular value. The curves have been plotted for two values of the coupling constant to the Higgs field, mf/mh ~ 1 and 
mf/mii ~ 4. Note the successive transitions between subsonic and supersonic behaviors according to the allowed jumps to 
the mass hyperbolae. However, the fermion vacuum mass limit does not allow the ultrarelativistic limit to take place in the 
perturbative sector, as it was the case for the lowest massive modes alone. In this case, the string dynamics follows a supersonic 
regime as soon as the first massive bound state is filled. 



25 



3. Discussion 



All these results have been derived without considering the back reaction effects induced by the trapped charge 



currents along the string [see Eq. (92)]. As was already discussed for the zero modes in Ref. ||25|], these currents 
yield back reacted gauge field, Bt and B^, which might modify the vortex background and the fermionic equations of 
motion. However, such perturbations of the Higgs and orthoradial gauge field profiles (see Fig. |]) can be neglected 
for i?fi?* and B^B^ small compared to the string forming gauge field BgB^ ^ ra^. Using Eqs. ( ^o|) and (|9^), the 
dimensionless charge currents associated with one massive bound state, and generating the dimensionless gauge fields 
Ba/nih, roughly read 

r - (104) 

with cj^ the fcrmion charges, i.e., cjFj^ or cj^-^. As a result, the back reaction on the vortex background is negligible as 
long as TO < 2iT^rj, which is clearly satisfied in the full perturbative sector. Moreover, since mh = Tj^fX, the previous 
derivations of the equation of state are also valid in the nonperturbative sector provided m < 27r^mh/\/A, and thus 
depend on the values of self-coupling constant of the Higgs field A, but also on the mass spectrum. As can be seen in 
Fig. 1^, the ratio m/TOf decreases with the fermion vacuum mass TOf , as a result m/TOh increases all the more slowly, 
which allows to have both TOf > TOh and to < l-K^rriY^j \f\. 

Moreover, in order for the new gauge fields Bt and B^ to not significa ntly modify the fermionic equations of motion. 



from Eq. (13), they have to verify qc^B^ <io ^ tojt. As a result, Eq. ( 104 ), and Ba/mh ~ j" yield the condition 



4- < 1- (105) 
r? C0 

As expected, it is essentially the same condition as the one previously derived for the zero modes alone ||2^. On 
the other hand, although the back reaction on the fermionic equations of motion can deeply modify the zero mode 
currents Q, since the massive bound states are no longer eigenstates of the ^^"f"^ operator, it is reasonable to assume 
that rather than modify their nature and existence significantly, the back reaction gauge fields may only modify their 
mass spectrum. In this sense, back reaction would indeed be a correction. 



VI. CONCLUSION 



The relevant characteristic features of Dirac fermions trapped in a cosmic string in the form of massive bound 
states have been study numerically, in the framework of the Witten model, and in the neutral limit. By means of a 
two-dimensional quantization of the associated spinor fields along the string world sheet, the energy per unit length 
and the tension of a cosmic string carrying any kind of fermionic current, massive or massless, have been computed, 
and found to involve as many state parameters as different trapped modes. However, in the zero-temperature limit, 
only two have been found to be relevant and they can be defined as the density numbers of the chiral zero mode 
exitations associated with the two fermions ^ and X coupled to the Higgs field. 

As a result, it was shown that the fixed trace equation of state no longer applies, as soon as massive states are filled, 
i.e., for energy scales larger than the lowest massive mode belonging to the mass spectrum. Moreover, the filling of 
massive states leads to a rapid transition from the subsonic regime, relevant with massless, or ultrarelativistic currents, 
to supersonic. Such properties could be relevant in vorton evolution since it has been shown that supersonic regimes 
generally lead to their classical instabilities [pl| . As a result, in the perturbative sectors for which TOf < TOh, the 
protovortons could be essentially produced at energy scales necessarily smaller than the lower mass of the spectrum, 
where the fermionic currents consist essentially in zero modes. In this way, vortons with fermionic currents could 
be included in the more general two energy scale models p^ . However, the present conclusions are restricted to 
parameter domains of the model where the back reaction can be neglected. Although it is reasonable to consider that 
the back reaction effects may simply modify the massive bound states through their mass values, their influence on 
zero modes are expected to be much more significant. In particular, the modified zero modes cannot be any longer 
eigenstates of the ^^^^ operator |Q, so one may conjecture that they acquire an effective mass, leading to massive 
states potentially instable for cosmic string loops 



26 



ACKNOWLEDGMENTS 



It is a pleasure to thank particularly P. Peter for many fruitful and helpfuU discussions. I also wish to acknowledge 
M. Lemoine, J. Martin, and O. Poujade for advice and various enlightening discussions. 



[1] T. W. B. Kibble, Phys. Rep. 67, 183 (1980). 

[2] D. P. Bennett, F. R. Bouchet, and A. Stebbins, Nature (London) 335, 410 (1988); D. P. Bennett and F. R. Bouchet, 
Phys. Rev. D 41, 2408 (1990); L. Perivolaropoulos, Astrophys. J. 451, 429 (1995); B. Allen, R. R. Caldwell, S. Dodelson, 
L. Knox, E. P. S. Shellard, and A. Stebbins, Phys. Rev. Lett. 79, 2624 (1997). 

[3] P. P. Avelino, E. P. S. Shel lard. J. H. P. Wu, and B. Allen, Phys. Rev. D 60, 023511 (1999); J. H. P. Wu, P. P. Avelino, 
E. P. S. Shellard, B. Allen, 



679 (1999). 



astro-ph/9812156 (1998); C. Contaldi, M. Hindmarsh, and J. Magueijo, Phys. Rev. Lett. 82, 



[4] F. R. Bouchet, P. Peter, A. Riazuelo, and M. Sakellariadou, Phys. Rev. Lett. , (to be published), |hep-ph/0005022| (2000). 
[5] Y. B. .Zel'dovich, Mon. Not. R. Astron. Soc. 192, 663 (1980). 
[6] A. Vilenkin, Phys. Rev. Lett. , 46, 1169 (1980); 46, 1496(E) (1980). 



[7] P. De Bernardis et ai, Nature (London) 404, 955 (2000); P. De Bernardis et at, |astro-ph/0105296 . 
[8] E. Witten, Nucl. Phys. B249, 557, (1985). 

[9] J. P. Ostriker, C. Thompson, and E. Witten, Phys. Lett. B 180, 231 (1986). 
[10] R. L. Davis, Phys. Rev. D 38, 3722 (1988). 

[11] R. Brandenberger, B. Carter, A. C. Davis, and M. Trodden, Phys. Rev. D 54, 6059 (1996). 

[12] R. L. Davis, Phys. Rev. D 36, 2267 (1987); C. T. HiU and L. M. Widrow, Phys. Lett. B 189, 17 (1987); M. Hindmarsh, 
ibid. 200, 429 (1988). 

[13] A. E. Everett, Phys. Rev. Lett. 61, 1807 (1988); J. Ambj0rn, N. K. Nielsen, and P. Olesen, Nucl. Phys. B310, 625 (1988). 
[14] B. Carter, Phys. Lett. B 224, 61 (1989); 228, 466 (1989); Nucl. Phys. B 412, 345 (1994); B. Carter, P. Peter, and 

A. Gangui, Phys. Rev. D 55, 4647 (1997); A. Gangui, P. Peter, and C. Boehm, ibid. 57, 2580 (1998). 
[15] P. Peter, Phys. Rev. D 45, 1091 (1992). 
[16] P. Peter, Phys. Rev. D 46, 3335 (1992). 

[17] A. Babul, T. Phan, and D. N. Spergel, Phys. Lett. B 202, 307 (1988). 
[18] J. Garriga and P. Peter, Class. Quantum Grav. 11, 1743 (1994). 

[19] P. Peter, Class. Quantum Grav. 11, 131 (1994); P. Peter and D. Puy Phys. Rev. D 48, 5546 (1993). 
[20] P. Peter, Phys. Rev. D 47, 3169 (1993). 

[21] B. Carter and X. Martin, Ann. Phys. (N.Y.) 227, 151 (1993); X. Martin, Phys. Rev. D 50, 7479 (1994). 
[22] X. Martin and P. Peter, Phys. Rev. D 51 4092 (1995). 
[23] B. Carter and P. Peter, Phys. Rev. D 52, R1744 (1995). 
[24] B. Carter and P. Peter, Phys. Lett. B 466, 41 (1999). 
[25] C. Ringeval, Phys. Rev. D 63, 063508 (2001). 

[26] S. C. Davis, W. B. Perkins and A. C. Davis, Phys. Rev. D 62, 043503 (2000). 

[27] N. K. Nielsen and P. Olesen, Nucl. Phys. B291, 829 (1987). 

[28] R. Jackiw and P. Rossi, Nucl. Phys. B 190, 681 (1981). 

[29] S. L. Adler and T. Piran, Rev. Mod. Phys. 56, 1 (1984). 

[30] S. A. Fulling, Phys. Rev. D 7, 2850 (1973). 

[31] B. S. Kay, Phys. Rev. D 20, 3052 (1979). 

[32] T. Goto, Prog. Theor. Ph ys. 46, 1560 (197 1); Y. Nambu, Phys. Rev. D 10, 4262 (1974). 
[33] P. Peter and C. Ringeval, hep-ph/0011308| (2000). 
[34] C. Ringeval, (in preparation). 



27 



